Reply to Comment on "Fluctuation-dissipation relations in the non-equilibrium 

critical dynamics of Ising models" 
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We have recently shown that in non-equilibrium spin systems at criticality the limit X°° of the 
fluctuation-dissipation ratio X(t, t w ) for t 3> tw 2> 1 can be measured using observables such 
as magnetization or energy [Phys. Rev. E 68, 016116 (2003)]. Pleimling argues in a Comment 
[cond-mat/0309652] on our paper that for such observables correlation and response functions are 
dominated by one-time quantities dependent only on t, and are therefore not suitable for a deter- 
mination of X°° . Using standard scaling forms of correlation and response functions, as used by 
Pleimling, we show that our data do have a genuine two-time dependence and allow X(t,t w ) and 
X°° to be measured, so that Pleimling's criticisms are easily refuted. We also compare with predic- 
tions from renormalization-group calculations, which are consistent with our numerical observation 
of a fluctuation-dissipation plot for the magnetization that is very close to a straight line. A key 
point remains that coherent observables make measurements of X°° easier than the traditionally 
used incoherent ones, producing fluctuation-dissipation plots whose slope is close to X°° over a 
much larger range. 
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In our recent paper [1] we analysed out-of-equilibrium 
fluctuation-dissipation (FD) relations in ferromagnetic 
spin systems quenched to criticality. One measures a 
response at time t to a perturbation at an earlier time 
i w , the waiting time, and compares with the correspond- 
ing two-time correlation. The FD ratio (FDR), X(t, t w ), 
then captures how much the response deviates from what 
would be expected in equilibrium. In many situations, 
T/X can be thought of as an effective temperature T e g 
governing the out-of-equilibrium dynamics [2]. In the 
context of critical dynamics the role of the long-time limit 
value X°° of X(t, i w ) for t ^> i w 3> 1 as a universal am- 
plitude ratio has been emphasized, see e.g. [3]. 

Exact results in [1] for the d = 1 case with T c = 0, 
showed that X°° is identical for all spin observables, 
whether one considers incoherent (short-range) observ- 
ables or coherent (long-range) ones such as the total mag- 
netization. For the coherent case, the limiting FD plot at 
long times is in fact a straight line, from which X°° can 
be determined trivially. Similar results were found for 
the case of bond observables, where the incoherent ob- 
servable corresponds to an indicator for a local domain 
wall while the coherent observable is the total energy. 
X°° was again found to be identical for these observ- 
ables. The advantage of using coherent observables is 
much more dramatic here: for the incoherent case the 
window in the FD plot where the slope is close to X°° 
shrinks to zero with increasing times, while for the co- 
herent observable one finds again a straight line FD plot 
of slope X°°. Numerical simulations in d = 2 strongly 
suggested that these results carry over to higher dimen- 
sions. In particular, FD plots for both the magnetization 



and the energy were numerically indistinguishable from 
straight lines. We also found the resulting values for X°° 
to be equal within numerical error, suggesting that there 
may be a well-defined effective temperature T c s for a 
broad range of observables. 

In a Comment [4], Pleimling argues that magnetiza- 
tion and energy are unsuitable for measuring X°° be- 
cause their correlation and response functions are dom- 
inated by one-time contributions depending only on the 
measurement time t. We show in this reply that Pleim- 
ling's remarks trivially apply in the regime t 3> t Wl but 
say nothing about the regime where t and t w are of the 
same order. It is in this regime that our numerical data 
were taken, and so they do carry non-trivial two-time 
information. Pleimling also argues that our results are 
in "marked contrast" with the results of renormalization 
group (RG) calculations [5] . We show explicitly that the 
RG results are in excellent agreement, predicting a limit- 
ing FD plot for the magnetization which is very close to 
a straight line. We also answer an "interesting question" 
raised, but left unanswered, by Pleimling. 

We begin by reviewing the construction of the FD 
plots from which we determine X(t,t w ), since Pleimling 
argues that the introduction of some one-time quanti- 
ties render our plots "unsuitable" . Consider a connected 
two-time correlation function C(t,t w ) = (A(t)B(t w )) — 
(A(t))(B(t w )}, with A, B two observables, and the conju- 
gate response R(t,t w ) — T S(A(t)) /ShB(t v/ )\ h =0 . Here 
fiB is the field thermodynamically conjugated to B and 
a factor of T has been included in the response. The 
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non-equilibrium FDR X(t,t w ) is defined via 

R(t,t w ) = X(t,t w )^-C(t,t w ). (1) 

This relation can be cast in terms of the step response 
x(t, £ w ) = f f dt' R(t, if), i.e. the response to a field hs 
switched on at t w and kept constant since: 

d d 
— x(t,tw) - -fl(Mw) = -X(i,t w ) — C(t,t w ) 

= X(t,t w )^-[C(t,t)-C(t,t w )}. (2) 

Two things are important to note. First, the correla- 
tion C(t, t w ) is a connected one, see the definition above. 
Second, it is physically sensible to compare the inte- 
grated response x(*>*w) = f t dt' R(t, t') to the integral 

fl dt' £rC(t, t') = C(t, t) - C(t, i w ) = AC(t, t w ), rather 
than just to C(t,t w ). These observations are irrelevant 
in the usual situation of incoherent spin observables, for 
which one-time correlations are constant, but are impor- 
tant in our case where they do change in time (a situation 
sometimes referred to as physical aging). This point is 
discussed in detail in Refs. [1, 6-8]. 

From (2) it follows that a parametric plot of x(t, t w ) 
versus AC(t,t w ) has slope X(t,t w ). This is obvious if 
i w is varied along the curve while t is held fixed. How- 
ever, if a series of such plots converges to a limit plot 
for t — > oo, then this limit plot and its slope X can 
clearly be obtained by varying either i w or t, as long 
as both times are large. For shorter times or if no 
limit plot exists, plots where t is varied and t w is fixed 
will not in general have a slope related to X, whether 
one plots x(t, t w ) versus C(t,i w ) or versus AC(t,i w ). 
However, in the simple case where X(t, t w ) is constant, 
one has from (2) that x(t, t w ) — XAC(f,t w ), so that 
a (AC, x)-plot does have the correct slope. A (C, x)- 
plot does not, on the other hand, since (dx/dt)(t,t w ) = 
X(d/dt)C(t,t) - X{dC/dt)(t,t w ) jt -X(dC/dt)(t,t w ). 
This lends further support to our choice of plotting x 
versus AC rather than C [9]. 

For systems where C(t,t) does not converge for 
t — > oo it is convenient to consider normalized 
functions x(Mw) = x(*> tw)/C(t, t) and C(t,t w ) = 
C(t,t w )/C(t,t) [1, 6, 8]. According to (2) these are also 
linked by 

J-X(t, t w ) = X(t, t w ) J-[l - C(t, tw)]. (3) 
ot w at w 

Again, X is the slope of a plot of \ vs 1 — C. If a limit 
plot is approached for large times, either t or t w can be 
varied to obtain this plot. Explicitly, if for large times X 
depends on t and f w only through C(t, t w ), the shape of 
the limit plot follows by integration of (3) as 

X(t,t w )= f dCX(C). (4) 

Jc(t,t v ) 
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FIG. 1: Normalized FD plot for magnetization in the d = 2 
Ising model at T c for times f w = 46, 193, and 720 (bottom to 
top). Curves have been vertically shifted by 0, 0.1, and 0.2 
for clarity. The convergence for large t w to an almost straight 
line of slope X°° ps 0.34 is evident. 



In equilibrium X(t, t w ) = 1 and one recovers the standard 
FDT relation x(t, t w ) = 1 - C(t, t w ). 

In [1], we showed FD plots for the total magnetiza- 
tion M = J2i s i (i-e. A = B — M above) and energy 
E = — Yluj) s i s j ior & d — 2 system of Ising spins Sj 
quenched to its critical temperature. These were pro- 
duced by varying t at several fixed t w , and without nor- 
malization. While a priori the slope of the plot does 
then not necessarily correspond to X, the numerical data 
for x(t,i w ) vs C(t, t) — C(t, t w ) fall on a straight line. 
Normalization only shrinks both axes of the plot in a t- 
dependent manner. The normalized plots will thus have 
the same slope, as shown explicitly in Fig. 1. The data 
clearly point towards the existence of a limit plot for large 
times which must be very close to a straight line. The 
asymptotic FDR A°°, which is obtained for t > t w > 1, 
i.e. C — > 0, is the slope at the end point of the limit plot 
(see the sketch in Fig. 2). Our data do not reach this 
end point, but RG calculations (see below) show that 
the slope should remain constant on approaching it. A°° 
can therefore be determined from the slope in the central 
part of the plot (i.e., the regime t > t w 3> 1). 

We now summarize the scaling relations used by Plcim- 
ling, taking as he did the case of the total magnetization 
as an example. For large times, one expects the two-time 
autocorrelation of M to scale as 

C(t,t w ) = tl +1 f c (t/t w ), (5) 

with a expressed in terms of standard critical exponents 
as a + 1 = (2 - rj)/z = (d - 2(3/v)/z. The scaling 
function fc decays as fc( r ) ~ t 8 f° r large r = t/t w 
with 0' w 0.19 for the d — 2 Ising model; in the limit 
r — » 1, fc(f) has to tend to a constant to have con- 
sistency with the scaling of the equal-time correlation 
C{t,t) ~ t (d-2p/y)/z = t a+x_ A simi i ar sca ling rela- 
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FIG. 2: Sketch of a limiting normalized FD plot (solid line). 
The asymptotic FDR X°° is the slope of the tangent at the 
top right end point of the plot (dot-dashed). This end point 
is at (1 — C, x) = (1, Y), where Y is the axis-ratio of the plot 
or, alternatively, the slope of the dashed line connecting the 
end point to the origin. 

tion holds for the response, R(t,t w ) — £^f R (t/t w ). As 
a result, X(t, t w ) becomes for large times a function of r 
only [3] ; this is confirmed by explicit RG calculations [5] . 
The normalized two-time correlation 

likewise only depends on r = t/i w only. Eliminating r, 
X can be expressed for large times as a function of C. As 
discussed above, it follows that a plot of x versus 1 — C 
must approach a limiting shape for large times; this is 
consistent with our numerical data in Fig. 1. Explicitly, 
if C(r) and X(r) are known then the limit plot is from (4) 

X(C) = ^ dr f - — J X{r) (7) 

where r(C) is the inverse function of C(r) and the minus 
arises because C(r) is a decreasing function. 

Pleimling deduces from (5) that, in the limit t — > 
oo at fixed t w , C(t,t) — C(t, t w ) « C(t,t) because 
C(t,t w ) <C C(t,t). This is correct, but not surprising. 
As stated above, one generically expects that t w ) — 
C(t,t w )/C(t,t) — > for t 3> t w . For the step response 
x(t,t w ) he shows similarly that this becomes indepen- 
dent of t w for t — > oo and grows with the same power 
law as C(t,t), so that the ratio x(£,i w ) = x(t, t w )/C(t, i) 
approaches a constant which we shall call Y [10]. Sum- 
marizing, for t — > oo at fixed t w , one has C — > and 
X — > Y. Referring to Fig. 2, Plcimling's statements thus 
fix a single point on the limiting normalized FD plot, 



FIG. 3: RG predictions for X(t,t w )/X°° as a function of 
the time ratio r = t/t„, for the two versions of the scaling 
functions F = 1 + e 2 AF (solid line) and F = exp(e 2 AF) 
(dashed line) . See text for details. 



namely its "end point" on the right. Geometrically, Y is 
the axis ratio of the FD plot. It is important to stress 
that Plcimling's reasoning says nothing about the rest 
of the limiting FD plot, which corresponds to the time 
regime where t and t w are of the same order: his limit 
t — > oo at fixed t w always implies the assumption t 3> tw- 
it is also clear from Fig. 2 that the axis ratio Y and the 
asymptotic slope X°° of the FD plot are not in general 
related. This answers the "interesting question" which 
Pleimling poses at the end of his Comment: Y and X°° 
are identical exactly if the limiting FD plot is a straight 
line. They have otherwise no reason to be related and 
there is no need to invoke RG calculations to "clarify" 
this issue. 

Pleimling's criticism would apply if, somewhat 
bizarrely, we had contrived only to collect data in the 
regime t > ( w . Such data would, in a normalized FD 
plot, fall very close to the plot's end point at (1 — C,x) — 
(l,Y). In an unnormalized plot, the t-dependent stretch- 
ing of the plot by C(t, t) would then indeed mean that 
the data trivially fall on a straight line. This line would 
be t w -independent and have slope Y rather than X°°. To 
check for such trivial behaviour, it is sufficient to normal- 
ize the data as explained above. We re-emphasize that, 
as Fig. 1 shows, our data arc not in the regime where 
such trivial behaviour is expected, covering a wide range 
of values of 1 — C and remaining well away from the end 
point of the plot. Our observation of a close-to-straight 
line FD plot is therefore not explained by scaling argu- 
ments, and remains highly nontrivial. This is transparent 
from Pleimling's own data [4] : one sees that his FD plots 
in Fig. 3 actually show data for which the t w -dependcncc 
of C(t,t) - C(t,t w ) (his Fig. 1) and x(Mw) (his Fig. 2) 
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FIG. 4: RG prediction for long-time limit of the normalized 
FD plot for the magnetization. Two versions are shown (solid 
and dashed) , corresponding to the two choices of scaling func- 
tion from Fig. 3, but are indistinguishable by eye. Both are 
extremely close to a straight line of slope X°° (dotted). To 
make the small differences visible, \ ~ X°°C is shown in the 
inset. 



is still significant. 

We now comment on Plcimling's statement that our 
numerical results arc in "marked contrast" to the RG 
calculations of [5]. These calculations give X(t, i w ) as a 
function of the time ratio r — t/t w in the form 



X(r) = X c 



F R (r) 
F dC {r) 



(8) 



where Fr and Fgc are appropriate scaling functions for 
J?(t, t w ) and (dC/dtv,)(t, t w ), consistent with the defini- 
tion (1). Both scaling functions are of the form F(r) = 
1+e A-F(r), within a second-order expansion in e = A—d. 
The extrapolation to d = 2, e = 2 therefore has a certain 
arbitrariness. To 0(e 2 ) in the RG calculation one could 
replace F(r) by e.g. exp(e 2 AF(r)). We show both ver- 
sions of the resulting RG predictions for X(r) in Fig. 3. It 
is clear that X(r) is close to X°° except for r w 1; where 
it does deviate, the RG predictions also become less reli- 
able. The near-constancy of X(r) already suggests that 
the FD plot will be almost straight. 

To find the limiting normalized FD plot predicted by 
RG explicitly, we combined the RG result for the scaling 
function Fgc(r) with the scaling exponents as quoted 
by Pleimling to obtain dC/dr and then used (7). The 
result is shown in Fig. 4 and demonstrates that the RG 
calculations predict a limiting FD plot which is extremely 
close to a straight line. Quantitatively, the plot is shifted 
upwards from a straight line of slope X°° by no more 
than Q.01X 00 ; its axis ratio Y therefore also lies no more 
than 1% above X°° . Contrary to Plcimling's remark, our 



numerical data are therefore entirely consistent with RG 
calculations. 

As a final point, we comment on our observation in [1] 
that the values of X°° are, to within numerical accuracy, 
identical for the magnetization (a spin observable) and 
the energy (a bond observable). While this may be sur- 
prising from the point of view of non-equilibrium critical 
dynamics [11], it is natural if one thinks of T c g — T c /X°° 
as an effective temperature which should govern the long- 
time non-equilibrium critical dynamics of all (or at least 
a broad range of) observables. This is supported by an 
analysis of the spherical model, where the values of X°° 
for spin and bond observables do indeed coincide. We 
will report on a more detailed investigation of this point 
in a future publication. 

In summary, Plcimling's criticisms of our method of 
measuring X°° using coherent observables do not ap- 
ply. His reasoning only addresses the limit^t 3> t w , 
where the normalized correlation function C(t,t w ) = 
C(t, t w )/C(t, t) is negligibly small, while our data are 
taken in a regime where C(t, t w ) is of order unity. This is 
most easily demonstrated using a normalized FD plot of 
X(t, t w )/C(t, t) versus 1 — C(t, t w ). Our observation that, 
for the magnetization in the d — 2 Ising model quenched 
to criticality, the normalized FD plot is close to a straight 
line therefore remains nontrivial, and is consistent with 
RG predictions. 

There are two key conclusions of our original study [1] 
which we have emphasized throughout this reply. First, 
FD plots for coherent observables are able to reveal non- 
trivial two-time dependences in non-equilibrium dynam- 
ics, and do so unambiguously when normalized. Sec- 
ond, FD plots for coherent observables typically have a 
wide range where their slope is close to the asymptotic 
value X°°. For measurements of X°° this makes them 
preferable to the traditionally used incoherent observ- 
ables, where this range shrinks to zero for long times. 

We are grateful to P. Calabrese for discussions. 
We acknowledge financial support from CNRS France, 
EPSRC Grants No. 00800822, GR/R83712/01 and 
GR/S54074/01, E.U. Marie Curie Grant No. HPMF- 
CT-2002-01927, the Glasstonc Fund, Nuffield Grant 
No. NAL/00361/G, Osterreichische Akademie der Wis- 
senschaftcn, and Worcester College Oxford. Some of the 
numerical results were obtained on Oswell at the Oxford 
Supercomputing Center, Oxford University. 



[1] P. Mayer, L. Berthier, J. P. Garrahan, and P. Sollich, 

Phys. Rev. E 68, 016116 (2003). 
[2] L. F. Cugliandolo, J. Kurchan, and L. Peliti, Phys. Rev. 

E 55, 3898 (1997). 
[3] C. Godreche and J. M. Luck, J. Phys. A 33, 9141 (2000). 
[4] M. Pleimling (2003), preprint cond-mat/0309652. 



■5 



[5] P. Calabrese and A. Gambassi, Phys. Rev. E 66, 066101 
(2002). 

[6] S. Fielding and P. Sollich, Phys. Rev. Lett. 88, 050603 
(2002). 

[7] A. Buhot and J. P. Garrahan, Phys. Rev. Lett. 88, 
225702 (2002). 

[8] P. Sollich, S. Fielding, and P. Mayer, J. Phys. Cond. 
Matt. 14, 1683 (2002). 

[9] When using normalized FD plots, the choice of AC vs C 
is immaterial, since AC(t,t v ) /C(t,t) = 1 — C(t,t v ) dif- 
fers from C(t,tw)/C(t,t) — C(t,t w ) only by a constant. 
[10] The limit t — ► oo here has to be understood as taken after 



the zero-field limit h — > 0. As Pleimling discusses, if one 
considers instead a fixed nonzero field then the response 
eventually becomes nonlinear, on a timescale for t which 
diverges for h — > 0. While correct, this observation is 
irrelevant for the present discussion: we have checked 
carefully that all our data are taken in the time window 
where the response is linear, as Pleimling's Fig. 2 also 
confirms. 

[11] RG calculations of X°° for the energy are in progress for 
the 0(n) model in order to clarify this point (P. Cal- 
abrese, personal communication). 



